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Recently, it has been conjectured that deconfining phase transition in SU(N C ) pure Yang-Mills 
theories is continuously connected to a quantum phase transition in softly broken M = 1 super 
Yang-Mills on 7? 1,2 x S 1 . We exploit this conjecture to study the strength of the transition and 
deconfining temperature as a function of the vacuum angle 8 in pure Yang-Mills. We find that the 
transition temperature is a decreasing function of 8 £ [0,7r), in an excellent agreement with recent 
lattice simulations. We also predict that the transition becomes stronger for the same range of 8, 
and comment on the 8 dependence in the large N c limit. More lattice studies are required to test 
our predictions. 



The presence of degenerate topological sectors in 
Yang-Mills theories entails adding a term Cq — 
—i 32 ^-i g -2 F^ iy F a,1L ' to the Lagrangian, where 9 is the vac- 
uum angle. In pure Yang-Mills or in the presence of 
massive fermions, both C and CV symmetries are lost. In 
QCD, the absence of CV violation puts stringent bound 
on the value of 9 < 10~ 10 . However, at finite temper- 
ature and close to the deconfinement transition, 9 can 
have spacetime dependence which gives rise to interest- 
ing physical phenomena in heavy ion collisions On 
the other hand, studying the nature of the deconfinement 
transition in Yang-Mills theories is plagued by strongly 
coupled physics. All this makes it very necessary to move 
to a simpler theory that is under complete analytic con- 
trol, yet resembles the original one. 

Recently, it has been argued that thermal phase tran- 
sitions in pure Yang-Mills SU(N C ) theories are continu- 
ously connected to quantum phase transitions in Yang- 
Mills theories formulated on R 1 ' 2 x S 1 , with one adjoint 
Weyl fermion of mass m and periodic boundary condi- 
tions [2, 3]. The mass term is crucial since the theory 
reduces to thermal pure gauge theory in the m — > oo 
limit, and to circle-compactified nonthermal M = 1 su- 
persymmetric theory as m — > 0. At small compactifica- 
tion radius L, the theory is weakly coupled and analyti- 
cally tractable. Near m = 0, the critical comactification 
radius L CT can be computed using reliable semi-classical 
methods. On the other hand, taking the limit m oo, 
one can interpret L CT as the inverse critical temperature 
T cr of pure Yang-Mills. 

At small L, the supersymmetric theory on R}' 2 x S 1 , 
with a gaugino mass term, dynamically breaks down (due 
to adjoint Higgssing in the S 1 direction) to [/(l)^" 1 at 
distances larger than 1/mw, where m\y — 2ir/(LN c ) is 
the W-boson mass. Then, the effective degrees of freedom 
are N c — 1 dual photons a (a N c — 1 dimensional vector) 
and <j) scalars which are the gauge field components along 
S 1 . We expand the fields <f> and a as <f> = j^-p+ j-b' , and 
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_ 2-nk+e -j-q-^ where g is the coupling constant, and 



the primed fields are the perturbations of <j> and a about 
their supersymmetric ground state. The ground state is 
labeled by an integer k = 0, N c — 1, the vacuum an- 
gle 9, and the Weyl vector p. For more details regarding 
the group theoretic definitions and normalization, see Q . 
The gauge holonomy (Wilson loop around S 1 ) in terms 

of the primed fields is D, = exp {i^-H ■ p + ij^H ■ b'J , 

where H are the Cartan generators H — (Hi, Hm c -i)- 
At b' — (supersymmetric vacuum), il obeys tr£7 = 0. 
We note that trfi is an order parameter of the Z^ c global 
center symmetry of the SU(N C ) theory which is unbroken 
in the confining phase. Thus, we have trfi = (trfi ^ 0) 
in the confining (deconfienment) phase. The total La- 
grangian governing the system is given by Eq. ([1]), with 
potential ([2]) which includes nonperturbative effects com- 
ing from monopole-instantons and bions. By studying 
this potential, we can determine the critical value of L at 
which the transition happens. 

Using the above mentioned conjecture, the authors in 
§ studied the phase transition of the pure SU(N C ) the- 
ories for all groups N c > 2 in the absence of the vac- 
uum angle 9. In this letter, we retain the vacuum angle, 
and study its influence on the critical temperature T cr as 
well as strength of phase transition. We find that T cr is 
a periodic function of 9 with period 2-7T, and is analytic 
everywhere except at (2n + l)7T,n G Z, as expected for 
physical observables [Ijj. For 9 in the range [0, tt), the crit- 
ical temperature decreases as 9 increases. This finding is 
in agreement with recent lattice studies that examined 
the dependence of the deconfinement temperature on 9 
On the other hand, the strength of the transition, 
quantified as the change of the trace of the Wilson line, 
increases with 9. For large N c , we find a weaker depen- 
dence of the deconfinement temperature and transition 
strength on the vacuum angle, due to N c suppression. 

The bion and monopole potentials and limit of 
metastability: The system effective Lagrangian follows 
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from the superpotential as given in [6], and describes 
the dynamics of SU(N C ) gauge group broken down into 
abelian [/(l)^ -1 subgroups upon compactifying on S 1 



C 



1 



m w 



12tt log(mvy/A) 



-Kp + Vr 



pert ■ 



(1) 



where A is the strong coupling scale, and V^ p and Vp er t 
are respectively the nonperturbative and perturbativc 
potentials. The perturbative part comes from Gross- 
Pisarski-Yaffe one-loop effective action for the holonomy 
Q. As was shown in Vp er t is suppressed by three 
powers of log(rnyp/A) compared to V np , and hence we 
neglect its contribution in our analysis. In the weak cou- 
pling limit, g 4 ^ c << 1, the nonperturbative potential is 
due to bions and monopole-instantons 0, and is given 
by 
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where ip = 2v £ +e , and the parameter k = 0, 1, .., N c — 1 
labels the vacuum branch, i.e. the branch with min- 
imum ground energy. We also introduced the dimen- 
sionless gaugino mass parameter c m = V^ on /V° ion . The 
bion and monopole coefficients V° ion and V^ on , expressed 
in terms of the physical mass mw and the strong scale 
A, are given by V° on = fl^log(^), and V£ on = 

9 ™A 3 i (mw_\ 
2tt m w 1U 6 V A J ■ 

First, we study the mass spectrum of the theory in 
the presence of the gaugino mass m and the 9 an- 
gle. To this end, we expand V np to quadratic order in 
b' and a'. We also rescale V and a' as {b' 2 ,a' 2 } — > 

^ lQ m" W ^ A ^ {^' 2 ' a ' 2 } ^° have a canonically normalized 
Lagrangian. Then, the mass term can be diagonal- 
ized by using the iV c -dimensional root vectors on — 
(0,0,..., 1 , -1,0,...) for i = 1,2,..JV C - 1 and 

i i+1 i+1 

o.n c = (— 1,0,0, This results in an extra massless 

particle that decouples from the rest of the spectrum. 
Performing these steps we obtain 
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and we notice that the 
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where m 2 , 

monopole potential causes a mixing between the b' and 



a' fields in the presence of a non- vanishing 6. To pro- 
ceed further, we use the Zn discrete Fourier transform 
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to find that the 



- - - i. o h 
quadratic part of the bion potential is diagonal in the 
{bp, dp} basis, while the quadratic part of the monopole 
potential has both a diagonal component oc cos ip as well 
as an off-diagonal component cx sin ip in the mass matrix. 
By performing an orthogonal transformation, we can go 
to new basis where the mass matrix is diagonal. The 
square of the mass spectrum in the new basis, as a func- 



tion of m, is rn± = 16m 2 



for 



^(^)±^sin 2 (^ 
p = 1, ...,N C — 1, and we have dropped the unphysical 
massless mode p = 0. We observe that the mass spec- 
trum is independent of the 8 angle. While this is true in 
the present case, restoring the g 2 N c term in the monopole 
potential in Eq. [2] will cause m\ to have an explicit de- 
pendence on 9, as was shown in the 577(2) case in 
Considering the values of m 2 _, we see that the center- 
symmetric vacuum becomes locally unstable for c m > c* n , 

where taking p = 1, we find c* n — 4 sin 2 (jj-j which is 
independent of 9. This corresponds to local instability at 
L* = A _1 a/ Itt-wt-t — — as one varies L for a fixed value 



of the guagino mass m. 

For any N c > 2 the non-perturbative potential con- 
tains cubic terms, and hence the transition is first order. 
In this the metastability limit above which 

the center-symmetric vacuum fails to be a local or global 
minimum. Unlike the metastable point c* m , the true crit- 
ical value depends on the value of the 9 angle, and is 
always smaller than c* m . 

Phase transition in SU{3): Before studying the 
phase transition, we must first determine the correct the- 
ory. Since for a given 9 we have k = 0,1,..., N c — 1 
branches, the correct theory is determined by select- 
ing the branch with the minimum vacuum energy. The 
theta dependence of the vacuum energy density is £ — 
min fc -Ysmm. cos (^ 27T ^ +e "j . For example, taking N c = 3 

the vacuum branch is k = for 9 G (— 7r, 7r), and k = 2 
for 9 G (tt,3tt). This implies that at 9 = n, the k = 
and k = 2 branches are degenerate 0, 0] (also, see 
for interesting discussions about the 9 term in a similar 
class of theories). 

Now, we turn to the study of the phase transition in 
SU{3). For 6 = 0, the V and a' fields do not mix to 
quadratic order. In this case, the expectation values of 
a' and b' 2 remains zero, and one can plot the potential V^p 
(a 2-dimensional plot) as a function of the deviation of 
the remaining holonomy from its center-symmetric value, 
i.e. b[, for different values of c m < c* n . Since c* m is the 
metastable limit, we can plot the potential for a range 
of c m < c* n = 4 sin 2 (§) =3 until we spot the true 
critical value c m at which the first-order transition takes 
place. Turning on a non-zero value for 9, will result in 
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FIG. 1: The N c = 3 and 8 = tt/2 effective potential V np 
(in units of V^ ion ) as a function of X'i f° r X2 = Xi — 0, and 
Xi=0. 



an admixture between a' and b 1 fields at the quadratic 
level, which does not enable us to see the phase transi- 
tion in 2-dimensional plots. To overcome this problem, 
we perform, as we mentioned above, an orthogonal trans- 
formation that brings us to new basis Xi,2,3,4 which are 
disentangled at the quadratic level. For 9 6 (— tt, 7t), 
the vacuum branch is k — 0, and the set of transfor- 
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Now, the 



is obtained by plotting V np against ip3 
for a range of c m < c* m = 3. This method is confirmed 
to a high precision for values of \9\ < 7r/2, by compar- 
ing it to the results of an optimization algorithm (errors 
less than 0.5%). Using this algorithm, we find that the 
expectation values of X2 an d X4 remain zero, while xi 
and X3 get expectation values such that |x3|/lxi| >> 1- 
As \0\ increases beyond tt/2 we find that |xi| gets expec- 
tation values comparable to those of \x3\1 and hence it 
becomes more difficult to see the first order transition on 
2-dimensional plots. Instead of trying to apply a second 
orthogonal transformation that disentangles xi an d X3 m 
the case of large values of \9\, we use the same optimiza- 
tion algorithm, that is checked for small values of \9\, to 
determine the true minimum of the potential. Multiple 
checks were carried out to confirm our results. 

The situation for SU(3) and 9 = tt/2 is shown in 
FIG. ([T]), and quite identical to 9 = case. Defining 
c** as the value of c m above which new metastable min- 
ima appears, we find that for c m < c** there is a unique 
center-symmetric global minimum (confining phase). For 
c*n < c m < Cmi there are three center- broken metastable 
local minima (only one is shown, the other two are ob- 
tained by applying Z3 rotations in the X3 — X4 plane), 
and a unique center-symmetric global minimum (confin- 
ing phase). For < c m < c* m — 3, the center-symmetric 
phase is a local minimum, and the three Z3 breaking min- 
ima are global degenerate minima. Finally, for c m > c* n , 
the center-symmetric point fails to be a local minimum. 
The critical temperature and strength of tran- 



FIG. 2: The critical normalized deconfinement temperature 
jcrfo) • ^ ne vacuum branch is k = for 8 G (— 7r,ir), and 
k = 2 for 9 £ (tt, 3tt) 



We find that LJ^- 



T ci - ( ) is a periodic 
function of 6 with period 2%, and cusps at 8 — (2n + 1)tt for 
any integer n, where two branches become degenerate. 
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FIG. 3: The normalized discontinuity of the holonomy 

r(0)/r(O) = A I j^L-trfi(0) I /A 1 3^trn(0)|. We take 8 in 

the range (— tt, 37r), and N c = 3. The vacuum branch is k = 
for 8 e (-7r,7r), and fe = 2 for e (7r,37r). We find that 
r(0)/r(O) is a periodic function of 8 with period 2n, and 
cusps at 8 = (2n + l)n for any integer n, where two branches 
become degenerate. 



sition in pure Yang-Mills: Now, we are in a posi- 
tion to determine the critical temperature for pure SU(3) 
Yang-Mills as a function of 9. To this end, we use 
the conjecture in 0, [3} that the thermal phase tran- 
sition in pure Yang-Mills theories is continuously con- 
nected to a quantum phase transition in softly broken 
W = 1 super Yang-Mills on R 1 ' 2 x S 1 . Using the relation 
c m = Con/K°ion = 167r 2 m/ (3A 3 L 2 Y c 2 ) , we see that at a 
fixed value of the gaugino mass to, the critical tempera- 
ture T CI = L- lcr is given by T cr (6>) = T oy /c%(6), where 
T = v /3A 3 Y c 2 /(16tt 2 to). In FIG. [2J we plot the normal- 
ized value of the critical temperature T" (9) / 'T CI (0) as a 
function of the vacuum angle 9 in the range (— tt, 27r). 
The critical temperature is a periodic function with pe- 
riod 2tt, and cusps at odd values of tt where two branches 
become degenerate and the CP symmetry is sponta- 
neously broken The normalized critical temperature 
decreases as 9 increases in the range (0, it). This is at- 



tributed to topological interference. The physics is that 
the effect of monopole-instantons is suppressed at 9 = tt 
with respect to 9 = by destructive interference among 
path histories. For small values of 9 (9 < 1), we fit- 



ted the ratio 
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K th = 0.003, and C - 



to the function 1 
= 2.09. 



K th 9 c , to find 



Such behavior is in excellent agreement with recent 
lattice studies that examined the 9 dependence of the 
deconfinement temperature in SU(3) pure Yang-Mills [B[ 
(also, see for a study of the 9 dependence across the 
deconfining transition, and also [l2[ for a review ). In the 
presence of 9, the Euclidean path integral representation 
of the partition function is not suitable for Monte-Carlo 
simulations because the measure is complex. The authors 
in [f| overcame this problem by studying the theory at 
imaginary 9, where the problem disappears, and then 
made use of analytic continuation to infer the dependence 
at real 9. By doing that, they were able to extract the 
9 dependence of T cr for small values of 9. They found 
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lattice 



= 0.0175. 



T cr (0) 

This is 5 times larger than the value of 1Z th . However, this 
should not come as a surprise since in general one does 
not expect an exact numerical match between a weakly 
coupled setup and pure Yang-Mills, a strongly coupled 
theory. 

The discontinuity at the first order phase transition 
can be quantified in terms of the change in the trace 
of the Wilson line. In the weak limit, g 2 N c /An « 1, 

we have trfi(0) = i (^) wj* (q H ■ &'(0)), where 



exp 
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N,. 



H-p 



is the center-symmetric holonomy. 
In FIG. [21 we plot the discontinuity of the absolute value 
of the holonomy, normalized to unity at 9 — 0, as a func- 
tion of 9 for N c = 3. We find that the change in the 
trace of the Wilson line has an opposite behavior to that 
of the critical temperature; it increases with increasing 
9. For small values of 0, we can fit A |trO(0)| /A |trfi(0)| 
to the function 1 + K9 C with K = 0.03, and C = 2. This 
behavior indicates that the strength of the phase transi- 
tion is enhanced for non-zero values of the vacuum angle. 
It is interesting to see if this trend can also be observed 
in lattice simulations by generalizing the study in [13j to 
non zero values of 9. 

We can also define the "latent heat" AQ of the sys- 
tem by identifying the quantity V np /L as the free en- 
ergy density F, and then employing the thermodynamic 
relation F = E — TS, where E is the internal energy 
density, T — l/L, and S — —dF/dT is the entropy 
density. Now, using T = Tq^/c{9), and realizing that 
the latent heat is the difference between the internal en- 
ergy at the symmetric and broken phases, we can show 
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0.4956 
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1.2441 


7r/12 


0.6735 


0.9107 


0.4953 


1.0680 


0.3790 


1.2447 


7T/6 


0.6724 


0.9114 


0.4946 


1.0683 


0.3786 


1.2462 


7r/4 


0.6706 


0.9128 


0.4935 


1.0689 


0.3778 


1.2486 


tt/3 


0.6680 


0.9144 


0.4918 


1.0692 


0.3767 


1.2517 


5tt/12 


0.6646 


0.9164 


0.4896 


1.0695 


0.3752 


1.2555 


tt/2 


0.6604 


0.9186 


0.4869 


1.0696 


0.3735 


1.2596 



TABLE I: The critical values of the gaugino mass parameter 



-trfi(e) 



c", and the Wilson line discontinuity T(9) = A 
for a range of 6 € [0, 7r/2], and iV c = 6, 7, and 8. 



(?' cr , 6' cr are the values of these quantities in the broken 
phase. We studied the quantity AQ(9)/AQ(0), to find 
that it has the exact same behavior of the trace of the 
Wilson line. 

The large N c limit and concluding remarks: For 

N c > 3, the nonperturbative potential V np is a function 
of N c — 1 holonomies. In this case, we use an optimization 
algorithm to find the minima of V np . In table (jT]) we 
list the critical values of the gaugino mass parameter c^j 
for a range of 9 6 [0, 7r/2], and for N c = 6,7, and 8. 
Similar to N c = 3 case, we observe the same pattern of 
the 9 dependence of the deconfinement temperature and 
the Wilson line discontinuity. However, as N c increases, 
the dependence of T cr and trfi on 9 becomes weaker. 
In the large N c limit, one expects that both quantities 
become independent of 9. Using the values of c^j listed in 
table H we can fit to the function 1 - K 2 9 2 /N V2 - 

U i 9 4 /N v ' 1 to find K 2 = 0.021, ft 4 = 0.0036, V 2 = 0.93, 
and T>4 = 2.36. However, we noticed that the values of 
T>2, and 2?4 get bigger as we restrict our fits to larger 
values of N c . It was numerically challenging to extend 
our analysis beyond iV c = 8 in order to check the stability 
of 2?2,4 as N c increases. 

In this letter, we have demonstrated the use of a class 
of Yang-Mills theories on R 1 ' 2 x S 1 as an analytic tool 
to learn about realistic gauge field theories. Interestingly 
enough, this class of theories is under complete control, 
yet resembles the original ones. It is the hope that, by 
contrasting the analytic predictions with the lattice data, 
these models will open the path into a new paradigm to 
tackle strong coupling problems. This work, along with 
0, B G3 , is a step toward exploring such a new path. 
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